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Abstract 



From the current cosmological observations of CMB and nuclear abundances 
we show, with an analytic procedure, that the total effective number of extra 
neutrino species AiV* ot < 0.3. We also describe the possible signatures of non 
standard effects that could be revealed in future CMB observations. This cos- 
mological information is then applied to neutrino mixing models. Taking into 
account the recent results from the SNO and SuperKamiokande experiments, 
disfavouring pure active to sterile neutrino oscillations, we show that all 4 
neutrino mixing models, both of 2+2 and 3+1 type, lead to a full thermaliza- 
tion of the sterile neutrino flavor. Moreover such a sterile neutrino production 
excludes the possibility of an electron neutrino asymmetry generation and we 
conclude that AiV* ot ~ 1, in disagreement with the cosmological bound. This 
result is valid under the assumption that the initial neutrino asymmetries 
are small. We suggest the existence of a second sterile neutrino flavor, with 
mixing properties such to generate a large electron neutrino asymmetry, as a 
possible way out. 
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I. INTRODUCTION 



Neutrino mixing is the simplest explanation of the data from the atmospheric [Q] and 
solar neutrino experiments 0|| , while alternative mechanisms are becoming more and more 
unlikely [|J. The results from the LSND experiment can be also explained by neutrino 
mixing ||. 

In this work we investigate the possible cosmological effects of neutrino mixing. The 
theory of big bang nucleosynthesis has been proposed since long time as a probe for particle 
physics ||[7]], but the systematic uncertainties in the measurements of nuclear abundances 
represent an obstacle for improvements. In the last two years different experiments confirmed 
the existence of acoustic peaks in the power spectrum of CMB temperature anisotropics |§ , 
from which it has been possible to measure, with improved precision, many different cos- 
mological parameters [p|-[T2]| . The precision of these measurements will be further improved 
by the new satellite experiments: the MAP satellite, already launched and on the way to 
its final orbit about the L2 Lagrangian point |IIJ and the PLANCK satellite whose launch 
is scheduled in the year 2007 [[14]]. These new observations represent a way to integrate the 
nuclear abundances observations partly overcoming the obstacle of systematic uncertainties 
and thus offering new opportunities to detect or constraint new physics in the early Universe, 
in particular the effects of neutrino mixing. 

In the section II we describe a simple new analytic and graphical procedure to confront 
a large class of possible non standard effects with the cosmological observations. We con- 
sider both the present situation, finding that the total effective number of neutrino species 
AiV* ot < 0.3, but we also point out which results from future observations could be inter- 
preted as signatures of non standard effects. In sections III-IV-IV we examine the specific 
predictions of those neutrino mixing models that have been proposed to explain the current 
data, including or excluding the LSND experiment and, with the new procedure, we confront 
them with the cosmological observations. We explain why the early Universe encounters dif- 
ficulties in detecting effects from three ordinary neutrino mixing models (section III), while 
it is emphasized the unique capability of the early Universe to probe a mixing with a new 
sterile neutrino sector even for very small mixing angles, otherwise out of the reach of Earth 
experiments (section IV). In the case of four neutrino mixing models (section V) we study 
the cosmological predictions using the results obtained in the simple active-sterile neutrino 
mixing and neglecting the possible presence of phases. We find the remarkable results that 
the new data from the SNO and SK experiments favour those four neutrino mixing models, 
both of '2+2' and '3+1 type', that lead to a final AiV* * ~ 1. In section VI we show how an 
additional sterile neutrino flavor could solve the disagreement with the cosmological bound if 
its mixing is able to generate a large electron neutrino asymmetry generation that produces 
a negative contribution to AA* ot . In section VII we conclude summarizing which are the 
possible signatures of neutrino mixing models that should be searched in future observation. 



II. COSMOLOGICAL OBSERVATIONS 
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A. Current constraints 



The recent observations of CMB anisotropics JT(]|-|r2[ provide a useful consistency test 
for the other cosmological observations. The interpretation of data depends on theoretical 
assumptions. Therefore, it is important that the simplest model used to fit the data, that 
makes use of 7 independent parameters, gives results that are consistent with the other 
observations. A combined analysis of the experiments allowing also for the presence of tensor 
fluctuations and a hot dark matter component by increasing the number of parameters to 11 



15| , does not show hints of the presence of such components and is remarkable that, when 
information from galaxy clustering is added, an upper bound of 4.4 eV on the sum of neutrino 
masses is found. Although these analysis support a cosmological consistency, one has to be 
aware that maybe we are excluding important possible effects that are still compatible 
with the data or that maybe some cosmological observations are affected by systematic 
uncertainties and are misleading us to wrong conclusions and to exclude important pieces 
of the picture. We will take the attitude to consider the simplest results as reasonable but 
at the same time we will check whether these assumptions are compatible with the neutrino 
mixing models that we will examine. 

In this section, we attempt to quantify the possibility that some non standard effects 
of BBN arise from neutrino mixing models. In order to test neutrino mixing models, 
CMB anisotropies are particularly important because they provide a measurement of the 
baryon content. This information has an important role in constraining the presence of new 
physics when is taken into account in models of Big Bang Nucleosynthesys (BBN). The 
BOOMERanG and DASI collaborations find an identical value [0)0] Q: 

(n b h 2 f MB = 0.022^3 (i) 



while the MAXIMA collaboration finds {fi b h 2 ) CMB = 0.0325lg;{£| at 95% c.l. pf. A 



combined analysis has been performed in ]T3| in which both hot dark matter and tensor 
fluctuations are allowed and it gives for the CMB alone at 95% c.l. (tt b h 2 ) CMB = 0.02±g;gf . 
If information from the IRAS PSCz survey on galaxy clustering is used then they find, at 68% 
c.l., (Vt b h 2 ) CMB = 0.0201^03 0- This result practically coincides with (|TJ) , even though 



different assumptions have been used. Therefore this seems quite a stable and reasonable 
value to be used for our analysis. 

The standard BBN model (SBBN) assumes the particle physics content of the standard 
model of particle physics (in particular zero masses and no mixing for neutrinos). Moreover 
it assumes that the neutrino distributions are described by the Fermi-Dirac ones with zero 
chemical potentials and with a temperature T v oc R^ 1 [0. The predicted primordial nuclear 
abundances are functions of the only parameter rj, the baryon to photon ratio, related to 
Vt b h 2 by the simple relation r] 10 = 10 10 rj ~ 273.6 Vt b h 2 . The value (111) for (fl 6 h 2 ) CMB 



1 We indicate la errors for all quantities unless differently explicitly indicated. More precisely the 
DASI experiment quotes at 68% c.l {Vt b h 2 ) CMB = 0.022^;^. 
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corresponds to have 0: 



V CMB = 6.0±Si (2) 
These predictions have to be compared with the measured values. A first group finds 'high' 



values for the primordial Helium abundance [18 



y p exp = 0.244 ± 0.002, (3) 
while a second group finds 'low' values |T9|]: 



y p exp = 0.234 ± 0.003. (4) 

At the moment there is a tendency to admit that there are systematic uncertainties in this 
kind of measurements and to unify the two ranges of values. We however prefer to continue 
to distinguish the two different measurements. 

The primordial Deuterium abundance is measured in Quasars absorption systems at high 
redshift. This kind of measurements gives the result f20f : 



(D/H) exp = (3.0 ± 0.4) x 10" 5 (5) 

We will not consider measurements of the primordial Lithium abundance since it is not fully 
understood whether we are really able to estimate how stellar processes could have modified 
it to the present. A test for the SBBN means to check whether the following conditions are 
fullfilled: 

Y SBBN {r] CMB )=Y e xp (g) 

(D/H) SBBN { V CMB ) = {D/H) cxp (7) 

The functions Yp BBN (rj) and (D/ H) SBBN (rj) do not have exact analytical expression but 
fits around r\ = 5 give the results [pl"|-|2"3"|j Q 

Y p BBN (v) - 0-2466 + 0.01 1h(t?/5) (8) 

(D/H) SBBN (ri) ~ 3.6 • 10~ 5 (r//5)- L6 (9) 
Using the CMB value (|2|) for rj, one finds that the SBBN predicts: 



2 From this moment we will always show values of rj in units of 10 10 , omitting the subscript '10' 
to simplify the notation. 

3 We are considering the neutrino heating from e + — e~ annihilations as a non standard effect (see 
discussion later on) and thus we are subtracting this contribution (AY p ~ 1.4 x 10~ 4 p4|] ) from 
the result found in Y^ BBN {r] = 5,r n = 887 sec) = 0.2467, where r n is the neutron life-time. 
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(10) 



SBBN 
'/highYp 


= 3.8 


f3.2 
-1.8 


SBBN 
'/lowYp 


= 1.4 


f2.1 
-0.8 


SBBN 
V(D/H) 


= 5.6 


f2.1 
-1.1 



(D/H) SBBN (r] CMB ) = (2.7 ± 0.7) x 10" 5 (11) 

If we compare these values with the experimental measurements we see that the SBBN is 
in agreement with the observations if high values of Y p are used, otherwise for low values of 
Y p there is a 4 a discrepancy. Such a comparison of SBBN with the observations can be also 
done saying that SBBN predicts, from the current Y p exp and (D/H) exp , the following values 
for r] (3a): 

(12) 
(13) 
(14) 

and comparing them with rj CMB the same previous conclusions follow. 

We want now to quantify the possibility that BBN is non standard and of course, in doing 
this, we will be particularly interested in those non standard BBN models that can result 
from neutrino mixing. In this case the possible non standard effects are of two kinds and 
quite well known. The first effect is the possibility that the number of energy density degrees 
of freedom g p = (30/tt 2 ) (p/T 4 ) differs from its SBBN value g s p BBN = (22/4) + (21/4) (T V /T) A 
before or during the BBN period. In this way the expansion rate and the standard BBN 
predictions for the primordial nuclear abundances would be modified 0. The change of g p 
can be expressed in terms of the (effective) extra number of neutrino species AN£: 

9 P = 9 S p BBN + l^ (I)' (15) 

From the definition of g p it follows that AN? is related to the neutrino energy densities by 
the following simple expression: 



A]V p = Px + Px _ 3 (16) 
Po 



where po — (77r 2 /120) T 4 and the 'X'-particles include the three ordinary neutrinos plus 
possible new species (we will be interested to possible new sterile neutrino flavors). Again 
we can make use of linear fits that account for the contribution of a non zero AN? in the 
BBN predictions for the primordial nuclear abundances Q: 



Y BBN ( V , AN p u ) ~ Y p SBBN (ri) + 0.0137 AN P V (17) 



4 The number 0.0137 can be inferred from the expansion given in [21] for 7] ~ 6. While the 
dependence of (D/H) on AN£ can be easily calculated considering that this abundance stays 



constant for r}j ' ^[g~ p = const [25 |. 
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(D/H) BBN (ri,AN£) ~ (D/H) SBBN (ri) (1 + 0.135 A^)°- 8 (18) 

A second class of deviations from SBBN are those related to distortions of electron neutrino 
and anti-neutrino SBBN distributions, given by the Fermi-Dirac distribution with zero chem- 
ical potential (the same for neutrinos and anti-neutrinos). In general deviations cannot be 
described in terms of a finite number set of parameters but by an infinite number of param- 
eters (the occupation numbers for each quantum state with a given momentum). However 
one can first calculate the change in Y p caused by these deviations and then normalize this 
change by introducing the quantity: 

AN U = \T P (ri,AN!,6f Ue>Pe )-Y p (ri,ANS)] 

0.0137 1 ' 

In this way one weighs the effect of distortions in terms of the presence of extra number 
neutrino species. A specific model of non standard BBN should be able to specify the 
deviations Sf Ue% p e {p,t) at each momentum and during all the period of BBN. However it 
has to be remarked that once that the neutron to proton ratio has frozen, the electron 
neutrino ditributions do not have any more a direct role in the nuclear reactions. Thus 
everything will depend only on the frozen value of n/p and still on AN£. This means that 
the Deuterium abundance will depend only indirectly on the electron neutrino distortions 
through the quantity ANp . Actually such a dependence is very weak and we will neglect 
it. Of course different 5f Ue ,u e can produce the same ANu and this degeneracy represents 
a lost of information The predictions of such non standard BBN models can again be 
compared with the experimental observations: 

Y p BBN ( v CMB , AN£, 8f UeiPe ) = Y^ (20) 

(D/H) BBN ( V CMB ,ANP) = (D/H) exp (21) 
The Y p measurement puts constraint on the quantity [J 



5 The only way to have more information on the 5f Ve ^ e would be to detect the electron relic 
neutrino distributions from which one could infer their values during BBN. Unfortunately relic 
neutrinos detection appears at the moment beyond the current observations but there are some 
interesting developments from the study of UHEz/ scattering on relic neutrinos and producing Z 
bursts [26|. 



6 Note there could be other kinds of non standard effects non considered here, like those ones 
associated with the possibility that during the BBN epoch there were baryon inhomogeneities on 



the scale of neutron diffusion length (see [27] and references therein). The quantity AN* ot — AN£ 
would assume a more general interpretation and depend also on other non standard parameters 
like the size of inhomogenenities and thus should be more generally indicated for example with 
paper we are interested to focus on non standard effects from neutrino mixing 
models and thus we completely neglect the possibility for these kinds of inhomogeneities but it is 
interesting that this procedure could be employed also in a different context. 
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AN tot = AN p + AN f Ve = Ay p SBJV /0.0137 (22) 
At 3a, assuming high values for F p exp , we find : 

AN tot = _o.3+o|, (23) 

while assuming low values we find: 

AA^ ot = —1.05 ± 0.75 (24) 

The Deuterium abundance provides a complementary information on AN£ and the compar- 
ison between the prediction and the observed value gives, conservatively at a 3a level, an 
upper bound on AN?: 

(AN P V ) BBN < 13 (25) 
while a lower bound is still not obtained with the current precision of measurements. The 



constraints (|23|), (E3) and are shown in figure 1, in a plot AA* ot — AN? 



B. Future observations and possible signatures 

It is interesting to see how one can expect that these constraints will improve from fu- 
ture CMB measurements of rj. The Planck experiment should be able to measure rj with a 



precision at the level of 1% or less p8|. In this way the uncertainties on the theoretical pre- 



dictions of the nuclear abundances, Y BBN and (D / H) BBN , will become negligible compared 
to the errors on the experimental values. 

Assuming that the future measured value rj CMB will correspond to the current central 
value of (D/H) in SBBN, rj = 5.6 (see Eq. flLl]) ), then the current Deuterium observations 
will constraint AN? to be < 4.0 (3 a) (the horizontal thick dashed line in figure 1) while 
still one does not get a lower bound []. The same exercise can be performed with Y p to see 
how the constraints on AA* ot would improve and the result is shown with vertical dashed 
lines in figure 1. This time the improvement is slight because the Y p abundance is much 
less sensitive to rj than (D/H) and the error on AA* * is dominated by the error on Y„ exp . 
The two constraints together, from (D/H) and Y p , give the gray region in figure 1 with 
thick dashed line contours. 

Differently it could happen that future CMB observations will indicate rj CMB > Tj^P^Lp 



see Eq. ([12])). In figure 1 it is shown, in light gray, the allowed cosmo logical region (at 



7 One finds AiV/? > —4 that is not particularly meaningful. At 2 a one gets ANp > —2, that 
implies the presence of at least one standard neutrino species. Note that a model in which a MeV r 
neutrino was decaying prior the onset of BBN was proposed to solve the BBN crisis from low Y p exp 



values [29|. In such a case one can get negative values of AN„ as low as —1. Nowaday such a model 
is very disfavored by the neutrino oscillations experiments but nevertheless it gives an example of 
why it is not meaningless to put a negative lower bound on AN{1, that moreover can be considered 
a sort of consistency check of the basic BBN assumptions. 
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3 a) in the plane AiV" — AiV£ for r] CMB = 7.0 (1% error) for a 'low+high' joint range of Y l 



v 

tot 



values. The SBBN would be ruled out at 3 a and negative AiV* ot < AN£ would be required. 
Therefore, in future, a 1% error measurement r] > 7 (or (Qt,h 2 ) CMB > 0.0256), will 
represent the opportunity to have a significant signature of non standard BBN effects with 
the current nuclear abundances observations. On the other hand the current allowed 3 a 
range of values of rj CMB , 3.6 — 9.9 (see Eq. (||)), excludes already now the possibility that a 
future 1% error measurement of r] CMB , with current Y p measurements, can give indications 
for positive values of AN£ Ve , that means AiV* * < AiV£, because it would require r] CMB < 2 
using high Y p values (see Eq. [12D and even lower values of r\ CMB using low Y p values. 

From the Eq. (0), one can see that from a 1% error measurement rj > 7.7 the 
Deuterium abundance would require also AN£ > (other than negative AN^). On the 
converse, for rj < 4.5, negative values of AN£ would be required. 

Another important improvement, from future observations of CMB anisotropies, will be 
the direct measurement of AN?. The presence of an extra radiative component changes the 
CMB spectrum, in particular leading to the enhancement of the height of first acoustic peak. 
At the present, a completely independent measurement of AN? from CMB anisotropies gives 
a very loose constraint (AN?) CMB < 19 (at 95% c.L, CMB alone) [0. However future MAP 
and Planck satellite experiments should reach a precision of 10 -3 — 10 _1 according to whether 
the information on the other cosmological parameters from other observations will be used or 
not and whether the CMB polarization will be measured or not . In figure 1 we indicated 
with horizontal thin dashed lines a realistic future constraint \(AN^) CMB \ < 0.1. It has to 
be said that the CMB observations will measure AN£ around the time of recombination 
and thus it could be in principle different from the value of AN£ during the earlier period 
of BBN if some intervening effect modified it. For example (AiV^) CMB can be higher of 
(AN£) BBN in the case of massive neutrino decays. In this case a comparison between the 
two quantities will test the 'relativity parameter' a oc m? u T [ |32]j . 

Their comparison could give also another result: (AiV^) CMB < (AiV^) BBN . This is 
possible only if one can say that (AiV^) BBN > 0. If one looks at the expressions (|3|) and 
(|14]), such a conclusion is possible if future 1% error observations will give r] CMB > 7.7. For 
example from r] CMB > 7.8 one can deduce (AN?) BBN > 0.2, while it can happen at the 
same time that CMB constraints AN^ MB < 0.1. This paradoxical situation could occur if 
AiV* ot is inhomogenous. We neglected a dependence of (D/ H) BBN on AiV", that means 
on Y p or, equivalently, on the frozen value of neutron to proton ratio. This because the Y p 
observations suggest that Y p cannot differ from the SBBN value so much to modify (D / H) 
in a sensible way, while the value of AN? is much less constrained and it can considerably 
alter the value of (D/H). However the observations measure Y p only within about 100 Mpc 
around us, while the (D/H) abundance is measured in the quasars absorption systems at 
much larger distances. Thus it cannot be excluded that Y p can be 'there' much larger than 
what we observe around us [B3J. The amplitude of CMB anisotropies exclude the possibility 



that this spatial variation can be due to a inhomogenous AN?, and thus it can be due only 
to a inhomogeneous AN„ that should be 'there' much larger (and positive) than around us. 
In this case the (D/H) nuclear abundance can be higher than what predicted by SBBN and 
compatible with r] CMB > r/ SBBN . Such a possibility should be however accompanied by the 
observation of dispersion in the (D/H) measurements in the range of values (1.8 — 3.6) x 10~ 5 . 
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Note that at the moment values of rj CMB > 7.7 are already excluded at 1.5<j and thus a small 
improvement in the measurement precision of r] CMB should be able to disfavour (or reveal !) 
such a situation. However only constraining the dispersion in the values of measured (D/H) 
can put more general limits on the presence of inhomogeneities in ANl Ve . 

Another important possibility is whether future observations will indicate AN£ > 0.3, 
because then, in order not to violate the bound AiV*°* < 0.3, one can conclude that there is 



a negative contribution AiV, 



C. Two special cases of non standard effects 



The SBBN corresponds, in the plane AiV* * — AN£, to the origin. One can consider 
the correction to the approximation of full neutrino decoupling at the time of electron- 
positron annihilations (implying T v oc R" 1 ). It has been shown that such a correction yields 
AiV* * ~ 0.012 and AN? ~ 0.034 p4|. Thus the predictions of nuclear abundances within 
the standard model of particle physics do not exactly coincide with those of SBBN. In the 
optimistic case that future CMB observations will be able to detect AN? as small as 0.01, 



then the small effect of neutrino heating should be distinguished [[31] and this would represent 
an important confirmation of the early Universe standard scenario below ~ lOMeV. 

A particular subclass of non standard BBN models, of the type considered here, is that 
one in which neutrinos have still thermal distributions but with non zero chemical potentials 
(fulfilling the chemical equilibrium condition /i Pq = —fi Ua )- 111 this particular case one has 
the following correspondence: 



30 /£, 

7 1 7T 



15 (§a 

7 { 7T 



AN/," 



■16 6 



(26) 



with £ a = /i a /T0. This kind of models have been largely studied in the literature since early 
times |1~7|J34]] . Constraints from the most recent cosmological observations have been also 
recently obtained in |]35 |. Their procedure put constraints on AN£, assuming that is equal 
in BBN and CMB epochs, and on 77 in a statistical combined way and taking into account a 
slight dependence of rj CMB on AN?. This allows to get more restrictive constraints but in a 
more specific context and at the expenses of physical insight. In our procedure we get more 
conservative constraints becasue of the poorest statistical procedure. On the other hand we 
gain more physical insight from an analytical procedure valid in a more general framework 
in which we distinguish AN? in BBN and CMB, the role of AiV> e is emphasized as we 
need for our purposes, and we find a bound on AiV*°* missing in [|35| . All these features are 
important for our following considerations. 



^The second relation is a good approximation for |£ e | <C 1 and AN„ <C 20. 
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III. THREE ORDINARY NEUTRINO MIXING 



With the exclusion of the LSND experiment, usually justified with the argument that 
it is the only experiment not yet confirmed by a second one, a three ordinary neutrino 
mixing can explain the current data from solar and atmospheric neutrino experiments. The 
three ordinary neutrino flavors eigenstates, \u a ) (a = e, //, r), are connected to three mass 
eigenstates \vi) with definite masses rrii (i = 1,2,3), by a 3 x 3 neutrino mixing matrix U: 

3 

K) = Y. U «M) (27) 

1=1 

The atmospheric neutrino experiments are then explained by the mixing of \i/ 2 ) and 
mass eigenstates with \m\ — m|| = <5m 2 tm ~ 2.5 x 10~ 3 eV 2 with a large mixing angle 
sin 2 2023 > 0.88 [IJ and with a negligible \U e3 \ <C 1 (as required by the CHOOZ experiment 



36fl ) that implies a small mixing angle #13. In this way the z/^'s are converted almost only 
to i/ r 's. The solar neutrinos are explained by the mixing of ui and z/2 eigenstates with 
\m\ — m 2 I = 5m 2 Q < 10~ 3 eV 2 . With the new data from the SNO experiment, large mixing 
angle solutions (sin 2 2^12 — 1) are also favoured |37| , |38|] . In this way the favoured three 
neutrino mixing models are those close to the the bimaximal mixing scenario [3S]. A three 
ordinary neutrino mixing does not have relevant effects on the cosmological picture and in 
particular on the quantities AiV£ and AiV/" e . It has been noted that a mixing of electron 
neutrinos with muon/tauon neutrinos during the period of freeze out of the neutron to 
proton ratio would exchange the abundances of the two types that are slightly different due 



to the different effect of neutrino heating [40]. In this way the effect of neutrino heating 



would change. However the neutrino heating effect is small like also the difference between 
the muon/tauon and the electron neutrino populations. Thus the bimaximal mixing would 
represent a correction of an already correcting effect. Therefore, at the moment, such a 
mixing model seems out of reach of cosmological investigation f\. 



IV. ACTIVE-STERILE NEUTRINO MIXING 

An explanation in terms of neutrino mixing of the solar and atmospheric neutrino exper- 
iments, together with the LSND experiment, implies three different scales of mass squared 
differences: Stjiq <C <5m 2 tm <C 5m 2 SND . This requires the existence of at least one new 
neutrino flavor [Ell] that has to be sterile in order to escape the constraint Ajf = 3.00 ± 0.06 



from the invisible decay width of the Z-boson ||42|| . A two neutrino mixing between one 
ordinary neutrino flavor v a and one sterile neutrino flavor v s is the simplest case of mixing 
involving new sterile neutrino flavors. With the new data from the atmospheric and the so- 
lar neutrino experiments, such a simple scheme seems to be excluded, as the three ordinary 



9 Cosmology is however useful to get information on the mass pattern, in particular structure 



formation and CMB observations are currently sensitive to a few eV's masses [15], while the 
PLANCK experiment will be be sensitive to a few O.leV's masses pSfl. 



10 



neutrino flavors appear to be mixed among themselves. However the solutions of the kinetic 
equations, necessary to calculate AN? and ANl ot , present many difficulties and this basic 
case represents an important starting point. Moreover it can represent a limit case for some 
of the possible sub-mixings within a realistic multiflavor mixing, as we will see in the next 
section. It is described by only two parameters, 5m 2 = m\ — m\ and s 2 = sin 2 2#o, where 
9q is the vacuum mixing angle. The \v%) mass eigenstate (l^)) corresponds to the ordinary 
(sterile) neutrino eigenstate in the limit of no mixing. The straightest cosmological effect is 
the sterile neutrino production with a consequent generation of a AN? that can be as high 
as 1 in the case of full thermalization. In doing these calculations one has to make some 
assumptions about the initial value of the ordinary neutrino asymmetries. We define the 
asymmetry of a lepton (baryon) particle species X as: 

Ny - Ny- 

- (28) 

7 

with Nx being the particle number per comoving volume and iV™ is the number of photons 
per comoving volume at an initial temperature T in ~ lOMeV ^> m e i/2 ~ 0.25 MeV. The 
(effective) total a-neutrino asymmetry is defined as: 

L<«) = L Ua + L Ve + L„„ + L Vr + Q a (29) 

with Q^t = -(1/2) B n and Q e = L e - (1/2) B n . For initial values |L (a) | < 
10 -6 (|5m 2 |/eV 2 ) 1 / 3 ('small' neutrino asymmetries) the effects on the oscillations can be 
neglected, while for much higher values an initial neutrino asymmetry can modify, usu- 
ally suppressing, the sterile neutrino production prior the onset of BBN For small 
neutrino asymmetries, the sterile neutrino production is given by |45|j4^| : 



N£ s = 1 - exp[-^(s 2 , |5m 2 |/eV 2 )] (30) 
The function g%{s 2 , |5m 2 |/eV 2 ) can be written in the form: 



gl(s 2 ,5m 2 /eV 2 ) = K a F«(s 2 ) s 2 yj \5m 2 \/eV 2 (31) 

where the subscript + (— ) stands for positive (negative) 5m 2 and K a ~ 657 (898) for a = 
e (/i, r). The function F±(s 2 ) is given by the following integral: 

t 2 

dt s 2 + d 2 t^ + (c±t^ (32) 

with c = cos 2^0 and d ~ 0.008(0.02) for a = e(/i, r). These results have been obtained 
within the static approximation |46| that neglects the MSW effect at the resonance. In 
the resonant case, for negative 5m 2 , this approximation holds only for very small mixing 
angles (s 2 10~ 4 ) and for small neutrino asymmetries |T7[] . Note that AN? = N£ s + 
(X^g ftT — 3), where the second term takes into account the depletion of ordinary neutrinos 
that is negligigible when the bulk of sterile neutrino production occurs before the neutrino 
chemical decoupling and in this case one has approximately AN? ~ Ng . This is verified 
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for \5m 2 \ > 10 4 eV that, for values of N Va > 0.01, corresponds to have small mixing angles 
s 2 < 10- 2 in the non resonant case and s 2 < 10 4 in the resonant case. In these regimes 
of small mixing angles, the function F±(s 2 ) is well approximated by its asymptotic value 
F±(0) and one has for a = e (fi, r): 

gl{s 2 < 10- 2 , l)/s 2 ~ K a F°(0) ~ 1.69 (2.33) 10 2 (33) 
g*(s 2 < 10-\ l)/s 2 ~ K Q (0) ~ 4.28 (2.27) 10 4 (34) 

In the non resonant case these analytical results agree very well with the numerical ones 
found in [^] Q for s 2 > 10~ 4 , while in |^9[ it is claimed that AN? is approximately three 
times lower. 

All these results have been obtained assuming small neutrino asymmetries. However in 
the resonant case, at small mixing angles, even if one starts with small neutrino asymmetries, 
a large a-neutrino asymmetry is generated around the critical temperature: 

T c ~ 15.0 (18.6)MeV (|5m 2 |/eV 2 ) 1/6 (2/y c ) 1/3 , (35) 



where y c is the critical (re-scaled) momentum [^,^,^T|]. The growth is first driven by the 
neutrino collisions that suppress the MSW effect. When the asymmetry reaches a value for 
which the interaction length at the resonance is larger than the resonance width, then the 



growth starts to be driven by the MSW effect p7| , |51| that can bring the z/ a -asymmetry up 
to a maximum value of 0.375. A remarkable feature is that the MSW effect is adiabatic 
for s 2 > 10- 9 (eVVl^m 2 !) 1 / 4 |If. Below this value the MSW effect b ecomes non adiabatic 



and ordinary neutrinos are not converted efficiently into sterile neutrinos anymore. However 
such a small value of the vacuum mixing angle represents by far the best example of how 
matter effects can enhance the mixing in vacuum, considering that in the Sun enterior the 
MSW effect occurs for s 2 > 10~ 4 loll. 



One would also like to know which is the upper limit on the vacuum mixing angle for the 
neutrino asymmetry to be generated. This is a point that has still not been investigated in 
the literature but we will see, in the next section, that it will prove to be very important for 
our considerations. Fortunately it is possible to get an analytic estimation. We will be par- 
ticularly interested to values of \5m 2 \ ~ 8m 2 LSND ~ 1 eV 2 . For these values one can use the 
Eq. ([34D to calculate the sterile neutrino production N£ . When the sterile neutrino produc- 
tion is negligible (N£ < 0.1) the value of the critical momentum is approximately given by 
the peak of Fermi-Dirac distribution: y c ~ 2 JE^,[HJ] . Once that the asymmetry generation 
has started the sterile neutrino production is suppressed in the collision dominated regime. 
This has the effect that the sterile neutrino production calculated by the Eq. ( |30"D is halved. 
Thus, taking into account this effect, one can easily calculate that the condition AN£ < 0.1 



10 In [Q| the results are presented for AiV", however for small enough mixing angles the contri- 
bution to AN£ Ue is negligible and a comparison is possible at least for the non resonant case. In 
the resonant case a comparison with the results of |4^| , at s 2 > 10 -4 , is not possible because these 
take into account also the negative contribution X^'^ N V0 — 3 to AN& and are performed in a 
quantum kinetic formalism that accounts also for the MSW effect. 
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corresponds to mixing angles s 2 < 0.52 (0.98) x 10~ 5 a/ eV 2 /\5m 2 \. When this condition is 
verified, together with the lower limit from the adiabaticity, the final value of the neutrino 
asymmetry is very close to the impassable limit corresponding to a situation in which, for 
an initial positive (negative) value of L^ a \ all anti-neutrinos (neutrinos) are converted into 
anti-sterile (sterile) neutrinos and thus |L;,J max = n Ua (n &a )/n 7 = 3/8 J50"| . Therefore the 
maximum value is also independent on the mixing angle in this range of values. When the 
sterile neutrino production becomes not negligible (AN? > 0.1), it has the effect to delay 
the asymmetry generation since the value of y c increases and therefore T c decreases. When 
y c becomes higher than ~ 10 the asymmetry generation at the critical temperature is driven 
by resonant neutrinos well in the tail of the distribution. Thus it is reasonable to think that 
when this happens the asymmetry generation mechanism will start to turn off. Unfortu- 
nately it is not easy to give an analytic description of this effect. However there is a much 
simpler reason for which the final value of the neutrino asymmetry has to decrease when 
the sterile neutrino production becomes not negligible. The reason is that the final value 
is reached during the MSW dominated regime that starts when the asymmetry has become 
large enough, during the collision dominated regime, that the neutrino and anti-neutrino 
resonances get completely separated and only anti-neutrino resonance can give a relevant 
effect, while the neutrino resonance is by far in the tail of the distribution if is initially 
positive, vice-versa if it is negative. In this way the MSW effect enhances the asymmetry 



to its maximum value | 50fl . However if sterile neutrinos have been produced during the 
previous collision dominated regime, not only ordinary anti-neutrinos will be converted into 
sterile anti-neutrinos, but also the already produced sterile anti-neutrinos will be converted 
back into ordinary anti-neutrinos. Thus the maximum value of the final neutrino asymmetry 
becomes 0: 

|£J max = ^(l-A^) (36) 

It will prove useful to assume, as upper limit on the mixing angle for the generation of neu- 
trino asymmetry, the value for which N£ > 0.9, corresponding to a final neutrino asymmetry 
at least one order of magnitude less than its maximum value 3/8. It is easy to calculate this 
value: 



s 2 ~ 0.5 (1) x 10- 4 ^eV 2 /\5m 2 \ (37) 

Let us discuss now the effects of a large neutrino asymmetry. The neutrino asymmetry 
generation yields two corrections to the N£ calculated, in the resonant case, from the Eq. 
(p0|). A first correction is due to the effect, just described, of suppression of the sterile 
neutrino production after the generation of the asymmetry and can be described by a cor- 
rective factor to the Eq. (RDf) that can be as low as 0.5, for N£ < 0.1 and that becomes 



quickly 1 (no suppression) for Ng > 0.1 [44|. A second effect takes into account the sterile 



n In H] it has been shown that the distribution function of sterile neutrinos produced during 
the collision dominated regime, for y c 3> 1, is given just by the equilibrium distribution times a 
coefficient a < 1 in a way that n Us /nl q = p Ua /pt q = Ny a = a. 
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neutrino production in the MSW dominated regime that results as an additive contribu- 
tion to N£ from the Eq. (^) that accounts only for the sterile neutrino production in the 
collision dominated regime. In the calculation of AN£ one has to take into account also 
the depletion of ordinary neutrinos. For —5m 2 <C 100 eV 2 the contribution to the sterile 
neutrino production from the MSW dominated regime occurs below the neutrino chemical 
decoupling and is compensated by an opposite ordinary neutrino depletion and thus there is 
no contribution to AN?. For higher values of —5m 2 ordinary neutrino are re-populated by 
the annihilations and this second contribution to AN£ can be as high as 0.4 |53[. However 
values of \5m 2 \ > 20 eV 2 are disfavoured by structure formation + CMB considerations (see 



for example ||15|| ). In any case the sum of the two contributions to the total AN£, from 



the two different regimes, cannot be much higher than 1. Thus the account of the neutrino 
asymmetry generation leads only to corrections to the calculation of AN£. 

In the case a = fi,r the contribution AN£ to ANjf 1 is the leading effect and we can 
approximately say that the accessible region in the plot AiV* ot — AN? lies along AiV* * = 
ANP for < ANP < 1 (see figure 2). 

In the case a = e, a large AN„ e can arise from two different processes. A first process is 
the u e , v e number density depletion that this time is a direct and relevant effect occurring for 
1 5m 2 1 < 10 _4 eV 2 and yields always a positive ANl"" that can be even higher than 1 0. The 
cosmological constraints are thus strongly strenghtened by the account of this effect |]56|j48[| 
and this can be seen in the plot AN^ ot — AN? considering that the region AiV* ot > AN? lies 
largely outside the cosmological allowed region (see figure 2). The second process is the 
generation of a large electron neutrino asymmetry in the resonant case and at small mixing 
angles. This time the sign is the same one of the initial that is observationally unknown 
and that could be predicted only within a full baryo-leptogenesys model and thus it can be 
both positive and negative. 

For negative values it is remarkable that the region AiV* * < AN? becomes accessible 
in the AiV* ot - AN? plot. In |k| it has been calculated that AiV"/" e can be as low as -1.4 
(for —5m 2 < 3eV 2 ) and in figure 2 one can see that for values AiV/" e < —0.3, there is 
compatibility with the region allowed by the low y exp values. Thus v e <-> u s oscillations 



provide a viable mechanism to solve the claimed SBBN crisis |57| , |50 



Another interesting possibility, shown in figure 2, is that v e <-> v s oscillations would 
be also able to reconcile possible future (1% error) values of -q CMB > 7 > ^sbbn ^ w ^ ^he 



12 There is a small positive contribution to AJV/" e at large mixing angles and |<5m 2 | < 10 4 due to 
a small depletion of v e number density induced by the much higher v a number density depletion 



[48|,|3|. 



13 For example in []54| it is shown that, for s 2 = 1 and 5m 2 ~ 3 x 10 8 eV 2 , the Y p production is 
0.02 higher than in SBBN, corresponding to AiV* ot ~ +1.5 and implying AN£ Ve at least as high 
as +0.5 (the value of AA^^ is not separately shown). Extrapolating to higher values of \5m 2 \ it 
seems also quite evident that much higher values of AJV" (3, 4, . . . ,?) and of AN„ e (2, 3, . . . ,?) 
are possible. This is confirmed by the results of a very recent work (5^] in which it is found that 
(A^ ot ) max ~ 6, implying AiV"/" e at least as high as ~ 5. 
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nuclear abundances observations 0. 

Still another interesting effect could be the possibility to generate the electron neutrino 
asymmetry in a inhomogeneous way. This would require the presence of small baryon 



number inhomogeneities |K| • This effect could produce inhomogeneous nuclear abundances 
that could have two kind of indications as discussed in the previous section: indirectly if one 
finds rf MB > r]f£§ N and (ANP) CMB < (AN„) BBN or, directly, if one finds a dispersion in 
the values of (D/H) measured from quasar absorption systems 0- 

Thus the generation of an electron neutrino asymmetry yields many interesting cosmo- 
logical effects but, within a two neutrino mixing scenario, it appears as a special possibility, 
considering that it requires a = e, negative 6m 2 and small mixing angles. However we saw 
that the generation takes place even for tiny values of the vacuum mixing angles and be- 
cause of this, the early Universe is the most sensitive way to probe small mixings with new 
sterile neutrino flavors. Moreover when considering realistic multiflavour mixing scenarios 
the conditions for the occurrence of an electron neutrino asymmetry generation can be more 
naturally satisfied. 



V. FOUR NEUTRINO MIXING 

Four neutrino mixing models represent the minimal way to explain, in terms of neutrino 
oscillations, all three anomalies including the results of the LSND experiment. These models 
are described by a 4 x 4 unitary mixing matrix U connecting the 4 mass eigenstates 
with definite masses mj, to the 4 flavor eigenstates \u a ) (a = e, /i, r, s): 

4 

K) = Y< U «M) (38) 

8=1 

There are different possible patterns for the mass spectra but all of them can be distin- 
guished in two types p^-p4]|. In a first type, the '3+1' models, the mass eigenvalue wi4 is 
separated by the other three, m*, by the LSND gap in a way that l^mfj ~ ^ m2 LSND- This 
case is a minimal modification of a three neutrino mixing model, since the introduction of a 
fourth mass eigenstate, to incorporate the LSND results, just perturbs the mixing among the 
other three explaining solar and atmospheric neutrino results. This means that the fourth 
eigenstate is almost coinciding with the sterile neutrino flavor Qk^l 2 — 1, \U as \ 2 Cl,«^s). 
In a second type of models, the '2+2' models, the spectrum splits in two nearly degenerate 
pairs with |<5m 2 | = 5m%., <5m 2 tm , separated by the much larger LSND scale Sm 2 LSND . In 



14 This possibility has been proposed in [pS[| , when the first data from BOOMERanG- MAXIMA 
were indicating rf MB = 9.0 ±1.4 Q. 

15 It is also interesting that inhomogeneous neutrino asymmetries, though on much smaller scales 
than those necessary to produce inhomogenous nuclear abundances, could be responsible for the 
generation of galactic magnetic fields and give rise to a detectable cosmological background of 
gravitational waves [61]. 
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this case the neutrino mixing matrix is very different from the case of three neutrino mixing 
models. 

There is an ongoing debate on which of the two types can better describe the experimen- 
tal data H63| , |65| - |6^1 . The new data from atmospheric neutrino experiments plus the inclusion 
of Tritium (3 decay data corner '3+1' models in two only allowed regions, at 99%c.L, around 
5m 2 LSND ~ 0.9 and 2eV 2 [67|. On the other hand the fact that both atmospheric |69[ 
Q and solar neutrino || data disfavour pure active-sterile oscillations, suggests, for '2+2' 
models, that the z/ e 's (for solar) and the z/^'s (for atmospheric) are converted into some 
admixture of both active and sterile neutrinos |7T|JBS[ |. Thus from solar and atmospheric 
neutrino experiments there is no evidence of the existence of sterile neutrinos and the sim- 
plest four neutrino mixing models, predicting pure active to sterile neutrino oscillations, are 
disfavoured. However there is not uncompatibility among the three experiments when the 
full range of possible four neutrino mixing models is considered. We will now study the 
cosmological effects of both '3+1' and '2+2' class of models. 

A. 3+1 models 

The '3+1' models can be distinguished in two classes, A and B, such that m 4 3> 
in A and m 4 -C in B. In the B case the three heavier mass eigenstates are almost 

degenerate with m 8 ~ sj Sm 2 LSND ~ 0.95 or 1.4 eV according to which of the two allowed 
islands is considered. The Heidelberg-Moscow experiment on (f3f3)o v decay puts restrictions 



on the B class ||72|| . We can make use of the results seen for v a <-> v s to get some simple 
estimations on the cosmological output of the two different classes. 

Class A. One has to consider the different possible ways of oscillations into the sterile 
neutrino flavor. The sterile neutrino flavor almost coincides with the fourth eigenstate but 
it is also slightly present in the other three eigenstates. The mixing between the three light 
eigenstates and the fourth eigenstate is set by Sm^ and, since it is positive, there is no 
neutrino asymmetry generation. The mixing of the three active neutrinos with the sterile 
neutrino can be described by three different mixing angles, sin 2 29 as ~ 4[/ 2 4 . For a = e, /i 
there are limits from the CDHS and BUGEY experiment for which sin 2 29 as < 10 _1 . For 
a = t we can assume the same limit. Thus from the mixing setted by Sm 2 LSND and using 
the Eq.'s fl3~0|) and ([33]) one can can see that there is a total thermalization (AN? = 1). The 



LSND experiment relates the two mixing angles in a way that sin 29 es x sin 29^ ~ 3 x 
10~ 4 . Therefore, even in the case of minimum sterile neutrino production, when sin 2 29 es = 
sin 2 29^ ~ 10~ 2 one has AN? ~ 0.9, very close to a complete thermalization. A mixing of 
the three active neutrinos with the sterile neutrino can be also driven by SrriQ and 5m 2 tm , 
since the sterile neutrino is also slightly present in the three light eigenstates. Now the sign of 
Sm 2 can be also negative and thus a neutrino asymmetry generation could occur in principle 
but the presence of a large sterile neutrino population, from the mixing set by Sm 2 LSND , will 
largely decrease the final value of the asymmetry, at least of one order of magnitude (see the 



16 c 



See also |f70|] for a critical discussion. 



16 



Eq. (j36D). In any case such a generation of neutrino asymmetry occurs for \5m 2 \ <C 10~ 2 eV 2 
and in this case the critical temperature would be lower than the freezing temperature of 
the neutron to proton ratio and would not affect BBN predictions in a way that | AN^ e | is 
negligible. Thus the only relevant cosmological effect is AN? ~ 0.9 

Class B. In this case the mixing of the three quasi-degenerate heavier eigenstates with the 
fourth eigenstate has a negative Sm^ ~ —$ m LSND- Therefore the sterile neutrino production 
is of resonant type and from the Eq. fl34|) with \Sm 2 \ ~ 1 eV 2 and sin 2 29 ~ 10 -2 — 10 _1 
one can see that again a complete thermalization would occur with AN? very close to 
1. In principle an electron asymmetry generation can also occur but the complete sterile 
neutrino thermalization has the effect to suppress completely the asymmetry generation 
mechanism and thus we can conclude that also in the B case \ANv\ <C 1 and therefore 

A7V tot = AN p ^ L 



B. 2+2 models 

These can be also distinguished in two classes, A and B. In the A (B) class the two lightest 
mass eigenstates, with masses m x and m 2 , explain solar (atmospheric) neutrino data while 
the two heavier, with masses and m^, explain the atmospheric (solar) neutrino data 
pB||>4] ]. Let us define simple limit cases in which the lightest and heaviest pair of mass 



eigenstates are made only of two flavor eigenstates, that means to consider a mixing matrix 
with two unmixed 2x2 blocks. Since the atmospheric neutrino experiments constraint 
the probability of — > v e conversions to be very small, then one has only four different 
possibilities: 

1. The rri3, mass eigenstates are made only of z/ M , v T and the mi, m<i mass eigenstates 
only of v e ,v s . 

2. The m 3 ,m 4 mass eigenstates are made only of v e) v s and the m^mj mass eigenstates 
only of v^v T . 

3. The m3, WI4 mass eigenstates are made only of u^, v s and the mi, m 2 mass eigenstates 
only of v ei v T . 

4. The m3,m4 mass eigenstates are made only of v ei v T and the mi,m2 mass eigenstates 
only of v s . 

Note that the 1 and 3 cases belong to the A class, while the 2 and 4 cases belong to the B 
class. A given neutrino flavor is only contained in one of the two pairs of mass eigenstates, 
that we call the normal couple, while it is absent in the other one, that we call the opposite 
couple. It is is simple to calculate the cosmological output since no neutrino asymmetry 
generation is possible (thus AN„ = 0) and AN? = 1 for the cases 3 and 4 and also for 
the cases 1 and 2 if the LMA solution is considered for the solar neutrinos otherwise 
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A LMA solution for a mixing v e 



is excluded by the Homestake experiment but it becomes 



possible is some mixture of v T is added to v s [38 or if Homestake is disregarded 
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AN? ~ 4 x 10" 4 for the SMA solution (sin 2 26 SMA ^ 10~ 3 , |5m 2 | SM A - 5x KT 6 eV 2 ||) 0. 

These simple four limit cases cannot explain the experiments for two reasons. The first 
reason is that in order to explain the LSND experiment the probability of — > v e conversions 
cannot vanish and the second reason is that the SK experiment and the SNO experiment 
H disfavour pure active to sterile oscillations. In order to explain the LSND experiment it 
is necessary that a small mixing between the heavy and light pairs of mass eigenstates is 
introduced in a way that there is a small contamination of v e and also in the respective 
opposite couple |63| , |64| . In order to explain the SK and SNO results one has to introduce a 



parameter that allows also for v e — > z/ MjT conversions [j in the case 1 and 2 and for z/ M — > v T 
in the cases 2 and 4. This is usually done introducing a mixing angle between the sterile 
and v^ T in a way that [y^ T , v s ) — > (v'u T , v' s ) = U(a){v^ T , v s ), where U(a) is a 2 x 2 rotation 
matrix |64||. In this way the 1 and 2 cases correspond to a = 0, while the 3 and 4 cases 
correspond to a = 7r/2 and for a = — > 7r/2 there is a continuous transformation bringing 
from 1 to 3 and from 2 to 4 . 

Class A. Let us first consider the transformation from 1 to 3. It is remarkable that 
when the condition for adiabaticity is satisfied for very small mixing angles sin 2 a > 
10 -9 (eV 2 /ISm 2 ]) 1 ^, a large muon-tauon neutrino asymmetry can be generated. This asym- 
metry generation can both suppress the sterile neutrino production from v e — > v s but 
also be partly transferred into an electron neutrino asymmetry yielding AN„ ~ —0.3 
or ANl"" ~ 0.1, according on the sign of the initial total asymmetry L^' T ^ ]73]| . However for 
sin 2 2a > 10~ 4 (see Eq. (|37|)), the mixing u^ T <-> u s with 5m 2 ~ — ^ m \sN di produces a ster- 
ile neutrino production AN£ > 0.9 that suppresses a large neutrino asymmetry generation 
and the final | AiV/"" | <C 0.1. The experimental data favour, like best fits, sin 2 a ~ 0.03 — 0.09 
and sin 2 a = 0.80 — 0.82 [^SJ Qand disfavour the possibility to have a (and also a' = ii/2 — a) 
smaller than 10~ 4 . Therefore the cosmological effects are a resonant sterile neutrino pro- 
duction with AN? ~ 1 and |AiV/" e | 1 for sin 2 a around the range 0.03 — 0.09 and a 
non resonant sterile neutrino production with AN? ~ 1 and |AA^/" e | = for sin 2 a around 
0.80 — 0.82. Note that the result is always AN? ~ 1, even assuming a SMA solution to solve 
the solar neutrino problem. In figure 3 the approximate accessible region for the class A 
'2+2' models and for sin 2 a > 10~ 9 (eV 2 /\5m 2 \f/ A is shown in the AA* ot - AN? plot with 
thick solid lines. The experimental results from the SNO and SK experiments corner it to 
the 'point' AN? = AiV* ot ~ 1, represented as a small circle. 

Class B. The other possibility is to consider a transition from the case 4 to the case 2 for 
a' = 7r/2 — a = — > n/2. Again when sin 2 2a' becomes > 10~ 9 (eV 2 / 1 5m 2 1 ) 1,/4 a neutrino 
asymmetry generation occurs and this inhibits the sterile neutrino production from z/ M <-> v s 



18 The SMA solution cannot give an electron neutrino asymmetry generation because 5m 2 is 
positive. 

19 With Vut we indicate a mixture of and v T . This further mixing has no relevance in cosmology, 
since the Vn and v T flavors cannot be distinguished. 

20 In the notation of 6£] coso = C23C24. 
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0. However again the SNO and SK experiments favour sin 2 2a' ^> 10~ 4 and the consequent 
large sterile neutrino production prevents a large neutrino asymmetry to be generated and 
again AN? ~ AiV* ot ~ 1. In figure 3 we again show a plausible accessible region in the 
AiV* ot — AN£ plot. Since in the case of asymmetry generation there are no calculations of 

£ 

\ANv e \, we show the most conservative region (between the dashed lines) assuming that 

£ 

can take all values between zero and the maximum possible value. This value 
corresponds to the case of an electron neutrino asymmetry generation in the limit of two 
neutrino mixing u e u s , for —5m 2 ~ leV 2 , as calculated in [53|. 



Let us try to quantify to which confidence level the results found in |68| constraint 
the possibility to have very small mixing angle sin 2 a < 10~ 4 or sin 2 2a' < 10~ 4 , that is 
equivalent to exclude the possibility of a neutrino asymmetry generation respectively in 
the A class and in the B class. This depends on which solution one assumes for the solar 
neutrino data. If one assumes a LMA solution then the best fit is for sin 2 a = 0.80 — 0.82 
or equivalently for sin 2 a' = 0.18 — 0.2. Very small values sin 2 2a' < 10~ 4 are excluded 
approximately at 95% c.l. If one assumes a SMA solution then the best fit is for sin 2 a = 
0.03 — 0.09. In this case very small values sin 2 a < 10 -4 are very slightly disfavoured and 
cannot be excluded to a significant statistical confidence level. However, from the reported 
values of Xmim the first case, assuming the LMA solution, is favoured compared to the 
second case, assuming the SMA solution, and thus values of sin 2 a < 10~ 4 are disfavoured 
approximately at the 90% c.l.. In the next year the SNO experiment should be able to 
constraint more significantly pure active to sterile neutrino oscillations in solar neutrinos 
and in particular the case when the SMA solution is considered, unless evidence for v e — > u s 
will be found. 

Therefore we arrive to the conclusion that current experiments favour those four neutrino 
mixing models, both of 3+1 and 2+2 type, in which the sterile neutrino flavor is brought to a 
complete, or almost complete, thermalization and no large electron asymmetry generation is 
possible in a way that the final result is always AN£ = AiV* ot ~ 1. Therefore, from the upper 
limit ANl ot < 0.3, current cosmological observations disfavour 4 neutrino mixing models. 
There are however some possibilities for which the cosmological bound could be evaded. 

• systematic uncertainties or statistical errors have been underestimated in Y p and/or 
rj CMB measurements. In the case of higher Y p and/or lower rj CMB then AiV* * = 1 
could be possible. For example one total extra neutrino species would be allowed at 
3 a in the case of underestimated systematic uncertainties if: Y p = 0.254 ± 0.002 and 
rj CMB unchanged as in the (0) or if rf MB = 3.5^Qg and high values Y p e * p are used; in 
the case of underestimated statistical errors if Y p exp = 0.244 ± 0.006 and rj CMB as in 
the (g) or if t] cmb = 6.0 ± 1.5 and high values of Y p cxp are used 0. 



21 Calculations of lAiV^I are missing in this case. 

22 These are qualitative estimations because we are calculating the 99% c.l. just multiplying by 
a factor 3 the error at 68%c.L, as for a gaussian distribution. This is a a very rough assumption 
when 5n/n is not <C 1 and a more elaborate statistical procedure should be used. 
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We assumed small initial neutrino asymmetries. If some unknown mechanism created 
large neutrino asymmetries (~ — 10~ 4 ) above T ~ 15MeV (the characteristic 
temperature for oscillations with \Sm 2 \ = 5m 2 LSND ~ 1 eV 2 ), then the sterile neutrino 
production would be completely suppressed [fl3[| . In this case the constraints that we 



obtained in section 2 should be applied to the values of AN? and AN„ associated to 
large neutrino asymmetries (see Eq. (^) ). 

We neglected completely the presence of phases in the four neutrino mixing matrix 
since we used the analogy with two neutrino mixing to calculate the cosmological out- 
put. The role of phases in cosmology has never been studied. A possibility could 
be that, when phases are taken into account, then the active-sterile neutrino mixing, 
even with large angles sin 2 2a > 10~ 4 , can generate a large neutrino asymmetry that 
suppresses the sterile neutrino production, or in the case of an electron neutrino asym- 
metry, yields a negative AN„ . This possibility should be verified in a full 4 neutrino 
mixing kinetic theory. 



VI. FIVE NEUTRINO MIXING 

If one assumes the existence of a mixing with a second light sterile neutrino flavor v s i 
then it is possible to evade the cosmological bound if the mixing generates a large neutrino 
asymmetry able to suppress the production of the first sterile neutrino and in the case of an 
electron neutrino asymmetry also to yield a large negative ANl" e . This new mass eigenstate 
should be added to the four neutrino mixing solutions that explain the experiments and that 
we described in the previous section. For convenience let us refer to the first neutrino flavor 
as the LSND neutrino. It is necessary to distinguish between models in which the LSND 
neutrino production is resonant and models in which it is not resonant. 

In the non resonant case, even though the asymmetry can start to be generated by the 
mixing with the s'-neutrino, afterwards it gets destroyed by the mixing with the LSND 



neutrino |43| , f44|1 . Thus the addition of a second sterile neutrino flavour to the A class '3+1' 
type models and to the '2+2' models that are close to the limit cases 2 and 3 cannot help 
to evade the cosmological bound. 

In the resonant case the generation of a neutrino asymmetry from the mixing with u s i 
is not obstacled by the mixing with the LSND neutrino. Thus any a-neutrino asymmetry 
generation has the effect to suppress the sterile neutrino production. However there cannot 
be a full suppression, because necessarily |<5m 2 s /| ~ |£m 2 J and the asymmetry generation 
starts when already about half of the sterile neutrino production occurred and the final result 
is AN? ~ 0.5. This is the only effect in the case of generation of a muon or tauon neutrino 
asymmetry and thus AiV* ot = AN? ~ 0.5. This means that adding a second sterile neutrino 
flavor to the '2+2' models 'close' to the limit case 1 (those for sin 2 a ~ 0.05) improves the 
agreement with cosmology but still not within 3a. In the case of an electron asymmetry 
generation one can have also a negative contribution from ANu e and the cosmological 
bound can be fully evaded. This means that the addition of a second sterile neutrino flavor 
makes possible to evade the cosmological bound only when it is added to four neutrino 
spectra of type '3+1' class B or '2+2' models close to the limit case 4 (sin 2 a ~ 0.80), 
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in which the LSND neutrino is present mainly in the light pair of mass eigenstates and 
an electron neutrino asymmetry can be generated. This possibility to evade the bound 
can be tested both with future f3f3 Qv decay experiments but also with future cosmological 
CMB observations that should find (AN?) CMB ~ 0.5. Moreover one should also have 
^cmb > yjSBBN ^ Du ^. considering the current error on D/H measurement, this possibility 

can be distinguished only if future CMB observations will give r\ CMB > 7.7 (at 3a). 



VII. DISCUSSION AND CONCLUSIONS 



We described an analytical and graphical procedure to search for non standard effects 
from nuclear abundances and CMB obervations. The recent measurement of the baryon 
content from CMB anisotropics improves remarkably the cosmological information on new 
physics. The present observations do not show evidences of the presence of non standard 
effects and constraints can be conveniently displayed in the AiV* * — AN? plot. 

However future measurements of 77 and AN? from CMB, together with the current mea- 
surements of primordial Helium-4 and Deuterium nuclear abundances, could provide some 
interesting signatures. Here is a summary list of the possible signatures as we found in the 
second section. 

1. If r] CMB > 7 then AiV* * < and a negative ANl Ve can be invoked. 

2. If r] CMB > 7.7 then also {AN£) BBN > 0. 

3. If (ANP) CMB > then (ANP) BBN > if one can exclude massive neutrino decays or 
other exotic effects intervening between the BBN and recombination epochs. 

4. If (AN?) BBN > 0.3 then, from the bound AA^ ot < 0.3, one can conclude that AA^" e < 
0. 

5. If the point 2 is verified but (AN£) CMB < (AN?) BBN , then this can be interpreted as 
a signature of inhomogenous AiV/"' 3 . This should be confirmed by inhomogeneities in 
(D/H) measurements that should be anyway searched independently on CMB obser- 
vations. 



We have applied these cosmological tools to the search of non standard effects from 
neutrino mixing. In the case of three ordinary neutrino mixing it seems not possible to 
find relevant cosmological effects. When a mixing with new light sterile neutrino flavors is 
considered, like the LSND experiment seems to suggest, then the early Universe becomes a 
powerful probe. We have shown how the SNO and SK experiments favour those four neutrino 
mixing models for which the sterile neutrino flavor is brought into thermal equilibrium or 
very close to it, implying that AN? ~ 1. At the same time a mechanism of electron neutrino 
asymmetry generation cannot be invoked to have a negative AN„ and thus in the end 
AA* * 1 fjne cosmological bound AAfJ ot < 0.3 is already quite conservative and future 
cosmological observations will be compatible with AiV" ~ 1 only if they will measure 
a value for r\ CMB that should be approximately half than the value measured by current 
observations or, alternatively, a value of Y p that should be about 0.01 higher. This of course 
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would mean that present observations are affected by large systematic uncertainties or that 
statistical errors have been largely underestimated. If one excludes such a possibility then a 
way out could be the presence of large initial neutrino asymmetries suppressing the sterile 
neutrino production. In this case the cosmological information can still be used to constraint 
the values of the neutrino asymmetries. Such a conclusion would have a quite remarkable 
impact on baryo-leptogenesis models. Another possibility is that phases in the neutrino 
mixing matrix could play an important role in the derivation of cosmological output and 
thus should be taken into account. Another intriguing possibility is to assume the existence 
of more than one sterile neutrino flavor. The new sterile neutrino flavor should be mixed with 
the electron neutrino flavor with the proper mixing parameters such that a relevant electron 
neutrino asymmetry is generated and both halves the sterile neutrino production and yields 
a negative AN„ e . This is possible only if the electron neutrino flavor is mainly present in 
the heavy mass eigenstates with m; ~ 1 eV. Therefore this scenario will be testable in future 
j3(3 0l/ decay experiments and with the cosmological tools that we described. 

This investigation thus shows that light sterile neutrinos in cosmology are now more 
constrained than before, because the possibility of a neutrino asymmetry generation in four 
neutrino mixing models is disfavoured within the statistical significance of the results from 
the SNO || and the SK [^] experiments. The result is that the sterile neutrino flavor, 
required by the LSND experiment, gets fully thermalized. Therefore the confirmation of 
the existence of light sterile neutrino flavors in next neutrino mixing experiments would be 
of great relevance for cosmology. Such a confirmation should come in next years by many 
planned experiments. In particular the MiniBooNE experiment should confirm or disprove 
the evidence of neutrino oscillations in the LSND experiment, while many other different 
experiments will be able to exclude exotic solutions to explain solar and atmospheric neutrino 
data. 
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Figure Captions 



Figure 1. Constraints on non standard BBN models from measurements of r]^ CMB \ Y p 
and (D/H). The solid vertical lines are the constraints (p3|) and (|24]) with the thick ones 
indicating the joint range coming from low+high Y p values. The horizontal solid lines is 
the constraint (|25|). The dark gray region is the allowed one by current observations. The 
dashed lines, contours of the gray region, are the constraints obtained neglecting the error on 
rj CMB in the BBN predictions and assuming a value -q( CMB ) = ^sbbn _ g ^ corresponding 
to (D/H) = 3.0 x in SBBN. The light gray region is the allowed region if one assumes 
rj CMB = 7 and low+high Y p exp range of values. The dotted horizontal lines are the realistic 
constraints that will be obtained on AN£ from future CMB anisotropy observations. The 
BBN from Standard Model of Particle Physics lies well within the circle centered around 
the origin. 

Figure 2. Accessible region for v a <-» v s in the plane AiV* 01 — AN£. The thick segment 
along the line AiV* ot = AN? corresponds to the case a = fi,r, while the striped regions 
are for the case a = e. The solid striped region is accessible in the case of u e , v e number 
density depletion below the neutrino chemical decoupling or in the case of negative electron 
asymmetry generation (AN„ e > 0) plus sterile neutrino production (AN£ > 0). The 
dashed striped region is accessible when a large positive electron neutrino asymmetry is 
generated. The thick dashed segment for AN? = corresponds to the region of mixing 
parameters for which the sterile neutrino production in the collision dominated regime is 
negligible. The possibility to have both an asymmetry generation and a sterile neutrino 
production (ANy ^ 0, AN? > 0) has not been studied in detail and there are only 
particular numerical examples. The dotted line is a simple interpolation between the two 
extreme cases AN? = and AN? = 1 that provides a reasonable approximation. 
Figure 3. Approximate accessible regions in the AiV* ot — AN? plot, in the case of four 
neutrino mixing models. The solid lines are for A class '2+2' models. The region between 
the dashed lines constraints B class '2+2' models. The dotted segment is for the A class 
'3+1' model. The small circle is for the B class '3+1' model and for all '2+2' models when 
the information from SNO and SK is used and it can be seen that it lies well outside the 
allowed cosmological region (in gray). 
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